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The  objective  of  this  test  was  to  record  the  frequencies  and  magnitudes  of  vibrJTIZ. 
and  shock  force  inputs  at  various  specified  locations  of  the  air  conditioner  during  the 
normally  expected  life  cycle  transportation  err/ironment. 


TEST  XTEM(B) 

Two  18,000  BIU  (horizontal,  carpact)  air  conditioners  (A/Cs)  were  mounted  in  a 
standard  S-280  shelter  in  a  manner  consistent  with  normal  (historical)  location  and 
hardware.  Upon  arrival  at  U.S.Army  Combat  Systems  Test  Activity,  Aberdeen  Proving  Ground, 
MD.  (USACSTA) ,  the  shelter  assembly  was  loaded  on  a  2  1/2  ton  truck  and  secured  using 
standard  tiedowns  and  procedures.  To  produce  a  worst  case  (maximum  shock  and  vibration) 
profile  for  truck  mounted  equipment,  the  shelter  was  enpty  except  for  the  installed  air 
conditioners.  The  truck  was  not  weighted  with  any  additional  ballast,  j 


TEST  PACTLITTES 

'  The  testing  was  conducted  on  road  courses  at  the  Munscn  And  Perryman  test  areas  of 
Aberdeen  Proving  Ground.  The  vibration  and  shock  environments  for  the  test  were  derived 
from  discussions  between  Belvoir  R  D  &  E  Crt  and  USACSTA.  The  various  courses  represent 
a  wide  selection  of  the  severety  of  possible  transport  environments  that  the  A/Cs  may  be 
subjected  to  during  a  normal  life  cycle.  The  limited  mileage  of  the  conducted  testing 
was  not  intended  to  represent  the  total  equivalent  wear  of  a  life  cycle  of  normal  use. 

Three  tri-axial  accelerometers  were  installed  on  the  roadside  A/C  at  the  locations 
marked  by  the  test  sponsor.  The  instrumented  testing  consisted  of  the  minimum  duration 
period  required  to  monitor  a  representative  sample  of  each  test  course  at  each  selected 
ground  speed.  Operational  speeds  of  the  test  vehicle  vere  controlled  by  a  pace  vehicle 
equipped  with  a  calibrated  fifth  wheel  driven  speedoneter.  Instrumented  testing  was 
ccnducted.cn  the  Munscn  area  improved  gravel  road,  Belgian  block,  two-inch  washboard, 
three-inch  spaced  bump,  and  six- inch  washboard  test  courses,  arrf  the  Perryman  PT1  level 
cross-ocxmitry  road  course.  Detailed  course  descriptions  are  contained  in  TECC24  TOP 
1-1-011. 

The  accelerometer  data  outputs  were  relayed  to  high  speed  magnetic  tape  at  an  onsite 
mobile  computerized  data  collection  van.  All  data  was  validated  at  the  tire  of  the  test. 
Data  reduction  to  the  format  requested  by  the  sponsor  was  accomplished  with  the  use  of 
special  ocrputer  programming  developed  at  USACSTA.  Details  of  the  instrumentation 
installation,  monitoring,  and  data  reduction  are  contained  in  the  engineering  report 
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REMARKS 


Prior  to  any  testing  the  operation  of  the  A/Cs  *sas  checked.  The  A/Cs  were  connected 
to  a  portable  electrical  generator  power  supply.  Operation  was  confirmed  by  listening  to 
the  fan  unit  operate,  and  confirming  cool  air  was  circulated  inside  the  shelter  and  heated 
exhaust  air  was  discharged  fraa  the  A/C  unit.  The  checks  were  performed  by  the  sponsor 
representative.  No  instrumented  measurements  were  required  for  the  operational  unit.  This 
brief  operational  capability  confirmation  check  was  performed  after  the  ccnpleticn  of  the 
instrumented  data  collection  for  each  different  test  course. 

In  addition  to  the  instrumented  mileage,  sate  rcninstrumented  mileage  was  accumulated 
to  complete  the  sponsor's  requested  testing.  The  operation  of  the  A/C  units  was  confirmed 
by  USACSTA  personnel  at  the  carpleticn  of  all  testing  in  the  same  manner  as  above.  At  the 
request  of  the  test  sponsor,  an  additional  150  miles  of  noninstrumented  transportability 
test  operation  was  conducted  cn  the  Perryman  PT1  cross-country  test  course. 

The  total  mileage  accumulated  during  the  testing  at  USACSTA  is  not  considered  a 
significant  portion  of  the  total  life  cycle  mileage  for  the  shelter  mounted  A/Cs  and  is 
not  considered  conclusive  of  any  findings  of  operational  capability  or  endurance. 

The  recorded  data  of  the  frequencies  and  magnitudes  of  vibration  and  shock  force 
Inputs  at  various  locations  of  the  air  conditioner  was  reduced  to  a  series  of  plots  of 
Rarer  Spectral  Density  (g^/Hz)  for  each  displacement  axis  at  each  mcnitored  location. 

Each  chart  is  identified  by  test  course  and  ground  speed  information.  The  data  is 
provided  to  the  sponsor  in  the  requested  formats.  The  data  may  be  used  to  produce  a 
random-cn-randcm  vibration  table  schedule  for  an  input  more  realistic  of  actual 
conditions  them  the  outdated  sine  sweep  schedule  currently  used  for  laboratory  test 
screening  and  development  of  A/C  units.  The  data  will  be  ccrparod  to  current 
specifications  to  determine  if  changes  in  future  production  requirements  should  be  made. 

A  detailed  presentation  of  the  test  setup,  instrumentation  and  data  interpretation  is 
presented  in  the  Engineering  Report  (End  1) . 


This  is  the  final  report  for  this  project. 
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SECTION  1.  EXECUTIVE  DIGEST 


1 . 1  SUMMARY 

The  shock  isolators  under  the  air  conditioners  did  rot  attenuate  the 
vibration  energy  to  the  air  conditioner  while  the  isolators  under  the  compressor 
did  reduce  the  energy  over  the  upper  parts  of  the  0  to  500  Hz  frequency  band. 

The  air  conditioner  compressor  was  subjected  to  shocks  caused  by  its 
impacting  other  components. 

The  vibration  test  schedules  developed  and  listed  in  Appendix  B  are 
dependent  upon  the  transport  scenario  given  in  ITOP  1-2-6C1,  1  March  1988.  This 
was  used  because  the  sponsor  was  unable  to  provide  an  item  specific  transport 
scenario.  Changes  in  the  total  transport  distance  and/or  the  percentage  of 
fatigue  inducing  road  surfaces  will  affect  the  laboratory  test  time. 

1.2  TEST  OBJECTIVES 

The  objectives  of  this  test  were: 

a.  To  determine  the  vibration  environment  during  road  operations. 

b.  To  develop  laboratory  vibration  test  schedules  to  simulate  the  field 
environment . 


c.  To  determine  the  ability  of  the  18,000  BTU  air  conditioners  to  operate 
after  exposure  to  the  field  vibration  environment. 


1.3 

TESTING  AUTHORITY 

None . 

1.4 

TEST  CONCEPT 

The 

The  18,000  BTU  air  conditioners 
testing  was  conducted  from  16 

were  subjected  to  road  vibration  testing, 
through  31  May  1990.  The  Vibration  Test 

Branch  of  the  U.S.  Army  Combat  Systems  Test  Activity  (USACSTA) ,  Aberdeen  Proving 
Ground  (APG) ,  MD  was  responsible  for  instrumenting  the  test  item,  conducting  the 
test,  and  acquiring  and  processing  all  test  data. 

1.5  SYSTEM  DESCRIPTION 

Two  50/60  Hz,  3-phase,  208 -volt,  18,000  BTU /hr  compact  horizontal  air 
conditioners,  KECO  Model  No.  F18H-3SB,  were  installed  in  the  front  of  an  empty 
S280(B)/G  shelter  (SN  138)  as  shown  in  Figure  B-l.  The  air  conditioners  weigh 
approximately  77  kg  (170  lb)  each.  A  rubber  grommet  was  located  above  and  below 
the  mounting  bracket  at  each  of  the  four  air  conditioner  mounting  bolts.  Four 
rubber  shock  isolators  were  installed  at  the  base  of  the  compressor. 
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2.1  ROAD  VIBRATION  TEST 


2.1.1  Objectives 

The  objectives  of  this  test  were: 

a.  To  determine  the  vibration  environment  during  road  operations. 

b.  To  develop  laboratory  vibration  test  schedules  to  simulate  the  field 
environment . 

c.  To  determine  the  ability  of  the  18,000  BTU  air  conditioners  to  operate 
after  exposure  to  the  field  vibration  environment. 

2.1.2  Criteria 

None . 

2.1.3  Test  Procedure 

The  empty  S-280  shelter,  housing  the  two  air  conditioners,  was  loaded  and 
secured  in  the  bed  of  an  M35A2C  2 -1/2 -ton  6x6  dropside  cargo  truck 
(VIN:  NK080J  0640-10781)  (fig.  B-l).  The  truck  was  equipped  with  size  9.00-20 
tires,  Load  Range  D,  inflated  to  345  kPa  (50  psi) .  The  roadside  air  conditioner 
(SN  873466)  and  mounting  hardware  were  instrumented  with  nine  uniaxial 
piezoresistive  accelerometers  arranged  in  three  orthogonal  axes  as  specified  by 
the  test  sponsor.  The  curbside  air  conditioner  (SN  873480)  was  not 
instrumented.  Photographs  of  the  accelerometer  locations  are  included  in 
Figures  B-2  through  B-4.  Table  B-l  contains  the  locations  and  specifications  of 
the  accelerometers  used  during  this  test. 

An  on-board  pulse  code  modulation  (PCM)  data  acquisition  system  was  used  to 
acquire  data  during  the  various  test  configurations.  A  detailed  description  of 
the  acquisition  system  is  presented  in  Appendix  G. 

Prior  to  testing,  an  electrical  calibration  was  performed  on  all  channels. 
The  accelerometers  were  calibrated  by  shunting  the  transducer  with  a  calibrated 
resistor  which  produced  a  differential  electrical  output  equivalent  to  a 
calculated  acceleration.  This  calibration  was  performed  in  both  the  positive 
and  negative  directions. 

A  linear  least  squares  curve  fit  was  performed  on  the  calibration  data  to 
determine  system  linearity  and  to  provide  a  quality  check  on  the  calibration 
data.  Any  channels  which  demonstrated  high  noise  levels,  DC  offsets,  and/or 
nonlinearities  could,  thereby,  be  detected  and  corrected  before  test  initiation. 
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2.1.3  (Cont'd) 

Vibration  data  were  acquired  while  the  air  conditioners  were  transported 
over  the  test  courses  at  the  speeds  listed  in  Table  2.1*1.  The  speeds  were 
determined  from  a  fifth  wheel  and  indicating  head  mounted  on  a  pace  vehicle.  An 
attempt  was  made  to  input  the  vehicle  speed  from  the  prime  mover's  speedometer 
into  the  data  stream.  This  was  not  done  because  the  vehicle's  speedometer  cable 
was  not  functioning  properly,  therefore  requiring  the  pace  vehicle.  Data  were 
recorded  for  approximately  31  seconds  during  each  different  speed  run  on  the 
test  courses.  The  short  length  of  two  of  the  test  courses  (Radial  Washboard  and 
Three-Inch  Spaced  Bump)  mandated  that  less  data  be  acquired. 


TABLE  2.1-1.  TEST  COURSES  AND  SPEEDS 


Speed 


Test  Course 

km /hr 

.  mph 

Belgian  Block 

2, 

16,  24, 

32 

5, 

10, 

15, 

20 

Two -Inch  Washboard 

16 

10 

Radial  Washboard 

24 

15 

Three -Inch  Spaced  Bump 

32 

20 

Six-Inch  Washboard 

4,  8 

2.5, 

5 

Perryman  Paved 

32, 

48,  64, 

80 

20, 

30, 

40, 

50 

Perryman  Cross-Country  No. 

1  8, 

16,  24, 

32 

5, 

10, 

15, 

20 

instrumented  air  conditioner 

was  operationally 

checked 

periodically  and 

at  the  conclusion  of  the  road  vibration  test. 

A  limited  amount  of  data  processing  was  performed  at  the  test  site  for  the 
purpose  of  data  verification  immediately  following  the  data  runs.  A  more 
complete  form  of  verification  and  all  data  analysis  were  performed  on  a  separate 
computer  system  after  the  data  acquisition  process  was  completed.  Data 
verification  and  analysis  procedures  are  presented  in  more  detail  in 
Appendixes  H  and  I . 

The  final  data  (after  verification)  were  processed  \jy  performing  the 
required  Fourier  transformations  and  creating  power  spectral  density  (PSD)  files 
for  the  individual  test  runs.  A  summary  of  the  parameters  used  to  compute  the 
PSD  functions  are  presented  in  Table  2.1-2. 


TABLE  2.1-2.  PSD  COMPUTATION  PARAMETERS 


_ Eaiameter _  Value 

Sample  rate,  Hz  2083.333 

Block  size  2048 

Frequency  bandwidth,  Hz  1.017 

Windowing  Hanning 

Number  of  linear  spectral  averages  32 
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2.1.3  (Cont'd) 

Amplitude  distributions  (rms,  +peak,  -peak,  +99.9  percent,  -99.9  percent, 
+99  percent,  and  -99  percent)  were  performed  on  all  channels  for  each  run.  The 
amplitude  distribution  data  were  compiled  (tables  3-2  through  B-7)  by 
histogramming  the  data  into  a  512 -bin  field  and  calculating  cumulative 
distributions .  The  average  value  for  each  channel  was  removed  to  account  for  DC 
offsets  in  the  instrumentation.  The  peak  values  presented  in  this  report  are 
the  true  peaks,  unlike  previous  Vibration  Test  Branch  reports  which  reported  the 
99.9-percent  value  as  the  peak  value  in  order  to  remove  extraneous  noise  from 
the  data.  The  current  data  verification  procedure  is  able  to  identify  noisy 
data  so  it  can  be  eliminated  from  the  analysis. 

The  laboratory  vibration  test  schedule  development  process  is  started  by 
grouping  (overlaying)  PSD  data  of  the  mounting  bracket  for  critical  speed  runs 
on  four  of  the  Munson  test  courses.  The  courses  and  speeds  selected  were  the 
ones  traditionally  used  in  the  development  procedure  and  are  listed  in 
Table  2.1-3.  These  PSDs  are  overlayed  to  provide  a  composite  spectrum. 

TABLE  2.1-3.  VIBRATION  SCHEDULE  TEST  COURSES  AND  SPEEDS 


Sneed 


Test  Course 

km /hr 

roph 

Belgian  Block 

32 

20 

Two- Inch  Washboard 

16 

10 

Radial  Washboard  ^ 

24 

15 

Three -Inch  Spaced  Bump 

32 

20 

The  overlay  process  was  accomplished  by  using  the  individual  average  plus 
one  standard  deviation  spectra  for  the  grouped  conditions.  The  use  of  the 
average  plus  one  standard  deviation  is  intended  to  add  some  level  of 
conservatism  to  the  data.  For  the  individual  spectra,  the  addition  of  a 
standard  deviation  has  the  effect  of  slightly  more  than  doubling  the  average 
value  while  the  addition  of  the  standard  deviation  during  the  overlay  process 
with  a  sample  size  as  small  as  two  runs  has  the  effect  of  enveloping  the  data. 

The  overlayed  data  were  then  exaggerated  to  reduce  the  laboratory  test  time 
to  a  more  reasonable  value.  A  real  world  exposure  time  for  cross-country 
operation  was  calculated  to  be  330  hours,  based  on  the  average  speed  of  the 
four  test  courses  utilized  and  the  expected  total  transport  distance  listed  in 
ITOP  1-2-601  dated  1  March  198S  40,000  km  (25,000  mi).  For  a  wheeled  vehicle, 
It  is  customary  practice  to  Ignore  the  vibration  levels  from  smooth  terrains 
since  those  levels  are  normally  lower  than  those  measured  while  the  vehicle 
operates  on  cross-country  surfaces.  Sixty-five  percent  of  the  total  distance  is 
assumed  to  be  on  cross-country  surfaces.  One-third  of  the  cross-country 
distance  is  assumed  to  be  of  the  severity  of  the  Munson  Test  courses.  This 
reduced  the  distance  to  8,700  km  (5,400  mi).  In  order  to  reduce  the  laboratory 
test  time  an  exaggeration  factor  of  2.0  was  applied  to  the  data.  This 
amplification  reduced  the  test  time  to  24.5  hours  for  each  of  the  three  axes. 
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This  is  still  s  relatively  long  time  for  a  laboratory  vibration  test;  however, 
the  exaggeration  factor  was  limited  to  2.0  because  it  is  unreasonable  to  amplify 
the  data  to  such  an  extent  that  the  test  level  will  exceed  the  yield  or  ultimate 
strength  of  the  material.  The  schedules  are  Included  in  Figures  B-5 
through  B-7.  Table  B-8  contains  the  breakpoints  for  these  laboratory  vibration 
schedules.  A  more  detailed  description  of  the  air  conditioner  environment 
schedule  development  is  presented  in  Appendix  J  . 

Frequency  response  (transfer)  functions  were  performed  by  dividing  response 
PSD  by  thair  corresponding  input  PSD  and  taking  the  square  root  of  the  result. 
This  is  discussed  further  in  Appendix  1. 

2.1.4  Test  Findings 

There  was  no  visible  damage  to  the  air  conditioners  as  a  result  of  the  road 
vibration  test.  The  air  conditioner  was  operational  throughout  and  at  the 
conclusion  of  the  road  vibration  test. 

The  acceleration  data  were  analyzed  in  the  form  of  amplitude  distributions 
(tables  B-2  through  B-7),  time  history  plots,  and  power  spectral  density  plots 
(fig.  B-8  through  B-l-*>). 

The  acceleration  levels  increased  on  each  course  as  the  vehicle  speed 
increased.  All  of  the  test  courses  with  the  exception  of  paved  road  had  a 
critical  speed  in  which  the  accelerometers  on  the  compressor  measured  shock 
loading.  A  typical  time  history  plot  showing  these  shock  pulses  is  included  in 
Figure  B-17.  These  shock  pulses  were  present  at  all  of  the  measured  speeds 
above  the  critical  speed. 

The  frequency  response  functions  indicate  that  the  vibration  energy  was 
amplified  through  the  shock  mounts  located  at  the  bottom  of  the  air  conditioner 
(fig.  B-18  through  B-20).  The  frequency  response  functions  between  the 
compressor  and  the  compressor  mounting  bracket  show  attenuation  in  the  vertical 
direction  throughout  almost  the  entire  0  to  500  Hz  band  (fig.  B-21) .  The 
transverse  and  longitudinal  transfer  function  plots  show  attenuation  above 
200  and  325  Hz,  respectively  (fig.  B-22  and  B-23). 

2.1.5  Technical  Assessment 

The  shock  isolators  under  the  air  conditioners  did  not  attenuate  the 
vibration  energy  to  the  air  conditioner  while  the  isolators  under  the  compressor 
did  reduce  the  energy  over  the  upper  parts  of  the  0  to  500  Hz  frequency  band. 

The  air  conditioner  compressor  was  subjected  to  shocks  caused  by  its 
impacting  other  components. 


2.1.5  (Cont'd) 

The  vibration  test  schedules  developed  and  listed  in  Appendix  B  are 
dependent  on  the  transport  scenario  given  in  HOP  1-2-601,  1  March  1988.  This 
was  used  because  the  sponsor  was  unable  to  provide  an  item  specific  transport 
scenario.  Changes  in  the  total  transport  distance  and/or  the  percencage  of 
fatigue  inducing  road  surfaces  will  affect  the  laboratory  test  time.  The 
control  accelerometer  for  the  laboratory  vibration  testing  should  be  placed  in 
the  same  location  as  where  the  data  were  obtained  (fig.  B-4) .  The  air 
conditioners  undergoing  the  laboratory  test  should  also  be  mounted  to  the 
vibration  exciter  in  an  identical  manner  as  those  during  this  road  test. 
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SECTION’  3.  APPENDIXES 
APPENDIX  B  -  TEST  DATA 


Table  No 

B-l 

B-2 

through 

B-7 

B-8 


_ Description _ 

Transducer  locations  and  specifications. 
Amplitude  distribution  data. 

Laboratory  vibration  test  schedule  breakpoints. 


El&VXg-No 

B-l 


B-2 

through 

B-4 

B-5 

through 

B-7 

B-8 

through 
B- 16 

B-17 

3-18 

through 

B-23 


_ Description _ 

Photograph  of  18,000  BTU  air  conditioners,  S280  shelter,  and 
M35A2C  truck. 

Photographs  of  accelerometer  locations. 

Laboratory  vibration  test  schedule  plots. 

Power  spectral  density  plots. 

Time  history  plot  showing  shock  loading. 

Sample  frequency  response  function  plots. 
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TABLE  B-l.  TRANSDUCER  LOCATIONS  AND  SPECIFICATIONS  18.000  BTU  AIR  CONDITIONER  TEST 
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TABLE  B- 3 .  ACCELERATION  AMPLITUDE  DATA 
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TABLE  B-5.  acceleration  amplitude  data 
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TABLE  B-7.  ACCELERATION  AMPLITUDE  DATA 
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TABLE  B-8 .  LABORATORY  VIBRATION  TEST  SCHEDULE  BREAKPOINTS 


u 


_ Vertical _ 

_ Transverse 

Loneitudinal 

Frequency 

Air.nl  itude 

Frequency 

Amplitude 

Freauency 

&rcp.Ul.yd‘ 

5.00 

0.5555 

5.00 

0.0555 

5.00 

0.1014 

6.00 

0.1526 

6.00 

0.0455 

6.00 

0.0811 

8.00 

0.2911 

7.00 

0.0479 

8.00 

0.2739 

9.00 

0.2161 

9.00 

0.1171 

9.00 

0.2672 

11.00 

0.0399 

10.00 

0.1200 

10.00 

0.3101 

13.00 

0.0409 

13.00 

0.0277 

14.00 

0.0558 

14.00 

0.0537 

14.00 

0.0298 

16.00 

0.0587 

15.00 

0.0419 

16.00 

0.0233 

17.00 

0.0771 

16.00 

0.0475 

19.00 

0.0402 

24.00 

0.0394 

21.00 

0.0243 

22.00 

0.0291 

27.00 

0.0425 

22.00 

0.0311 

22.00 

0.0206 

28.00 

0.0315 

26.00 

0.0268 

23.00 

0.0245 

30.00 

0.0481 

28.00 

0.0163 

25.00 

0.0142 

32.00 

0.0357 

30.00 

0.0249 

27.00 

0.0182 

35.00 

0.0681 

34.00 

0.0319 

28.00 

0.0153 

41.00 

0.0698 

35.00 

0.0419 

29.00 

0.0191 

47.00 

0.0425 

37.00 

0.0379 

34.00 

0.0142 

49.00 

0.0648 

39.00 

0.0463 

41.00 

0.0245 

56.00 

0.0234 

40.00 

0.0379 

45.00 

0.0149 

57.00 

0.0366 

42.00 

0.0475 

49.00 

0.0338 

69.00 

0.0300 

46.00 

0.0335 

54.00 

0.0110 

70.00 

0.0385 

49.00 

0.0475 

59.00 

0.0264 

74.00 

0.0265 

51.00 

0.0268 

63.00 

0.0186 

81.00 

0.0447 

54.00 

0.0327 

68.00 

0.0444 

122.00 

0.0058 

56.00 

0.0163 

69,00 

0.0277 

130.00 

0.0096 

58.00 

0.0268 

81.00 

0.0613 

150.00 

0.0101 

62.00 

0.0176 

85.00 

0.0284 

161.00 

0.0061 

64.00 

0.0261 

91.00 

0.0383 

178.00 

0.C023 

67.00 

0.0163 

149.00 

0.0047 

199.00 

0.0089 

70.00 

0.0282 

155.00 

0.0074 

215.00 

0.0014 

71.00 

0.0167 

195.00 

0.0022 

226.00 

0.0047 

81.00 

0.0255 

217.00 

0.0024 

270.00 

0.0024 

102.00 

0.0189 

224,00 

0.0069 

287.00 

0.0032 

109.00 

0.0104 

240.00 

0.0052 

310.00 

0.0012 

128.00 

C . 0296 

265.00 

0.0093 

368.00 

0.0135 

135.00 

0.0148 

314.00 

0.0037 

457.00 

0.0C05 

149  00 

0.0275 

320.00 

0.0069 

481.00 

0.0038 

168.00 

0.0107 

346.00 

0  0034 

500.00 

0.0018 

182.00 

0.0204 

371.00 

0.'X>61 

193.00 

0.0199 

386.00 

0.00.8 

226.00 

0.0013 

402.00 

0.0058 

260.00 

0.0159 

448.00 

0.0006 

500.00 

0.0026 

500.00 

0.0012 
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Figure  B-2.  Accelerometer  location  at 


B-3.  Accelerometer  location  at  compressor  top. 


Accelerometer  location  at  air  conditioner  mounting  bracket. 


18K  ITU  AIR  CONDITIONER  VII  SCHEDULE  VERT 


Figure  B-5.  Vertical  laboratory  vibration  test  schedule. 
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Figure  B-6 .  Transverse  laboratory  vibration  test  schedule. 


Figure  B-7.  Longitudinal  laboratory  vibration  test  schedule 
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Figure  B-8.  Power  spectral  density  plot,  Belgian  block,  32  km/hr. 


RUH  006  (T)  COMPRESSOR  10TT0H  (flVE) 


Figure  B-9 .  Power  spectral  density  plot,  Belgian  block,  32  kra/hr. 
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Figure  B-ll.  Power  spectral  density  plot,  Belgian  block,  32  km/hr. 


density  plot,  Belgian  block,  32  km/hr. 


Figure  8-13.  Power  spectral  density  plot,  Belgian  block,  32  kin/hr. 
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Figure  •>- 15 .  Power  spectral  density  plot,  Belgian  block,  32  km/hr. 
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Tirae  history  plot  showing  shock  loading. 
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Figure  B-19.  Frequency  response  function  plot,  Belgian  block,  24  kra/hr. 
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Figure  B-21.  Frequency  response  function  plot,  Belgian  block,  24  km/hr. 
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Figure  B-23 .  Frequency  response  function  plot,  Belgian  block,  24  km/hr. 


APPENDIX  G.  PULSE  CODE  MODULATION  DATA  ACQUISITION  SYSTEM 

Each  of  the  transducers  were  connected  to  a  USACSTA  signal  conditioner 
which  provided  excitation,  amplification,  calibration,  and  low-pass  filtering. 
The  accelerometers  were  low-pass  filtered  at  500  Hz.  Interchangeable  six-pole 
Butterworth  filter  pads  were  installed  in  the  signal  conditioner  with  the  filter 
frequency  being  the  nominal  frequency  at  which  the  transducer  response  is 
attenuated  by  3  dB.  A  block  diagram  of  the  PCM  data  acquisition  system  is 
included  in  Figure  G-l. 

The  output  from  each  channel  of  the  signal  conditioner  was  connected  to  an 
EMR  Model  372-03  Pulse  Code  Modulation  (PCM)  Encoder.  The  encoder  multiplexed 
the  incoming  analog  signal  which  was  then  sampled  at  2083.333  samples  per 
second  per  channel  by  a  sample- and-hold  amplifier,  digitized  by  a  10-bit 
successive  approximation  analog-to-digital  converter,  and  converted  to 
nonretum-to-zero-level  (NRZ-L)  code  for  transmission.  A  programmable  read-only 
memory  was  used  for  the  encoder  format  control.  The  encoded  PCM  data  were  the 
inputs  to  a  Conic  Model  C7L  510  230.9  MHz  Transmitter  for  transmission  to  the 
remote  data  handling  facility.  The  signal  conditioner,  encoder  and  transmitter 
were  mounted  in  the  cab  of  the  prime  mover.  The  transmitted  NRZ-L  code  was 
received  by  a  Scientific  Atlanta  Series  420S  Receiver  and  recorded  on  a 
Honeywell  Model  1-1  PCM  Tape  Recorder  along  with  voice  annotation  of  the 
individual  test  runs  and  IRIG-B  time  code.  Simultaneously,  the  PCM  pulse  train 
was  passed  to  a  Loral  Instrumentation  ADS-100  PCM  Decommutator. 

The  ADS-100  system  consists  of  sevjral  modules.  The  input  buffer  and  bit 
synchronizer  modules  recovered  the  serial  pulse  train  from  the  data  link  noise 
and  disturbances.  The  frame  synchronizer  and  data  distributor  modules 
demultiplexed  the  serial  pulse  train  into  16-bit  words.  A  parallel  input  module 
was  used  to  input  digital  IRIG-B  time  code  into  the  ADS-100.  The  compressor  and 
direct  memory  access  (DMA)  modules  performed  data  sorting  and  transfer  to  the 
host  computer,  an  HP2?  -E  series  minicomputer.  The  data  were  temporarily 
stored  on  the  HP  system  uisk  and  later  transferred  to  digital  tape  to  provide  a 
permanent  storage  medium.  The  ADS-100  system  was  independent  of  host  control; 
however,  software  residing  on  the  minicomputer  controlled  the  hand  shaking 
between  the  ADS-100  and  the  HP21MX-E  during  data  acquisition. 
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Figure  C-l.  Block  diagram  of  PCM  data  acquisition  system. 
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DATA  VERIFICATION 


The  data  verification  that  was  performed  consisted  of  searching  the 
individual  channels  for  frame  errors,  wild  points  and  DC  shifts.  A  frame  error 
occurs  whenever  the  PCM  stream  is  interrupted,  which  induces  artificial 
high  -level  spikes  into  the  data  stream.  Data  approaching  full  scale  values,  in 
addition  to  frame  errors,  are  flagged  by  the  frame  error  check  program,  because 
the  program  defines  a  frame  error  as  being  any  data  value  which  is  greater  than 
98  percent  of  full  scale.  The  same  program  checks  the  data  for  wild  points  and 
DC  shifts.  A  wild  point  is  just  like  a  frame  error  except  its  amplitude  is 
smaller  and  does  not  approach  full  scale.  The  program  defines  a  wild  point  as  a 
positive  or  negative  change  of  more  than  100  resolution  steps  from  the  previous 
point.  A  DC  shift  occurs  whenever  the  data  average  changes  rather  abruptly  and 
then  remains  constant  after  the  change.  The  program  defines  a  DC  shift  as  a 
shift  of  more  than  25  resolution  steps  from  the  previous  average,  with  the 
average  being  based  on  approximately  1/6  second.  These  errors  are  obvious  when 
plotted  in  the  form  of  a  time  history,  but  the  quantity  of  data  collected  is 
such  that  viewing  it  all  as  a  time  history  is  impossible.  For  this  reason, 
the  frame  error  check  program  is  utilized,  and  the  locations  of  all  such  errors 
on  the  system  disc  are  printed  out  for  a  permanent  record. 

As  a  further  step  in  the  data  verification  process,  acceleration  amplitude 
distribution  data  were  compiled  by  histogramming  the  data  into  a  512-bin  field 
and  calculating  cumulative  distributions .  Table  H-l  shows  sample  acceleration 
amplitude  distribution  data.  The  average  value  for  each  channel  was  removed  to 
account  for  DC  offsets  in  the  instrumentation.  The  percentile  columns  represent 
the  percentage  of  time  the  data  was  in  that  particular  range.  For  example, 
99.9  percent  of  the  time,  the  data  value  for  channel  1  was  less  than  1.02  but 
greater  chan  -0.71.  The  units  for  the  accelerometer  amplitude  distribution  data 
are  g's. 

TABLE  H-l.  SAMPLE  AMPLITUDE  DISTRIBUTION  DATA 


Run  Number,  Test  Course,  Speed,  etc. 


Description 

nas 

♦Peak 

-Peak 

+99,?* 

-99.9% 

+99% 

-??* 

+?9* 

-99% 

Channe 1 

1 

0.28 

1.25 

-1.05 

1.02 

-0.71 

0.68 

-0.55 

0.37 

-0.36 

Channel 

2 

0.31 

1.32 

-1.02 

1.08 

-0.78 

0.71 

-0.62 

0.39 

-0.42 

Channe 1 

3 

0.23 

1.34 

-0.83 

0.96 

-0.65 

0.66 

-0.50 

0.29 

-0.27 

Channel 

4 

0.22 

1.22 

-0.82 

0.94 

-0.58 

0.62 

-0.42 

0.26 

-0.26 

Channel 

5 

0.24 

1.04 

-1.02 

0.85 

-0.64 

0.59 

-0.49 

0.33 

-0.31 

Channel 

6 

0.36 

1.50 

•1.57 

1.19 

-0.92 

0.81 

-0.73 

0.46 

-0.46 

Channel 

7 

0.39 

1.66 

•1.37 

1.30 

-0.97 

0.86 

-0.81 

0.50 

-0.49 

Channel 

8 

0.21 

1.27 

-0.83 

0.95 

-0.60 

0.64 

-0.44 

0.24 

-0.24 

H-l 


The  program  which  performs  the  amplitude  distribution  analysis  and  creates 
tables  such  as  the  one  above  performs  a  number  of  data  validity  tests  for  each 
channel  and  provides  messages  such  as: 

a.  Channel  inactive. 

b.  Data  one-sided. 

c.  Data  noisy. 

d.  Data  clipped. 

e.  Large  rms  value. 

f.  Large  DC  offset. 

g.  No  data  spread. 

h.  Shock  present  in  data. 

If  a  large  number  of  samples  are  processed  and  if  the  dynamic  range  of  the 
data  i.i  small,  a  bin  overflow  can  occur.  This  happens  when  more  than 

32,767  samples  occur  in  any  one  bin  and  usually  appears  as  an  inactive  channel 
(rms  set  to  0.0).  Although  the  amplitude  distribution  program  is  not  foolproof, 
and  the  rules  which  determine  data  validity  are  somewhat  arbitrary,  the  program 
provides  a  very  useful  tool  for  "quick  look"  verification  and  analysis. 

The  test  data  were  then  checked  for  stationarity  or  time  invariance.  Data 
are  collected  and  analyzed  for  short  periods  of  time  wherein  seasonal 

differences  in  performance  and  long-term  vehicle  degradation  have  no  influence. 
Hence,  a  single  time  history  record  should  adequately  define  the  distribution  of 
the  data.  Multiple  sample  records  are  used  in  all  steady- state  analyses  under 

the  assumption  that  the  data  are  stationary.  This  provides  a  broader,  more 

dependable  sample  population  and  statistical  base  for  analyses  such  as 
acceleration  amplitude  distributions  and  acceleration  power  spectral  densities. 
The  test  for  stationarity  is  based  on  the  following  assumptions: 

a.  The  data  to  be  analyzed  are  not  of  transient  nature  (e.g.,  single 
mechanical  shock) . 

b.  The  sample  record  taken  will  reveal  the  non- stationary  character  of  the 
random  process  in  question. 

c.  The  sample  record  is  very  long  relative  to  the  lowest  frequency 
component  in  the  data. 

d.  Any  nonstationarity  of  significance  will  be  revealed  by  changes  in  the 
rms  value  of  the  data  with  time. 
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Because  the  objective  is  to  measure  short  term  rapidly  changing  phenomena 
rather  than  seasonal  differences  or  vehicle  degradation,  it  is  assumed  that  a, 
b,  and  c,  for  this  purpose,  are  true.  A  nonparametric  statistical  test,  called 
a  run  test,  is  used  to  check  assumption  d.  The  run  test  consists  of  computing 
the  rms  value  for  a  series  of  subrecords  from  the  time  history  (e.g.,  30  rms 

values  over  1-second  intervals  for  a  30-second  record),  finding  the  median  value 
of  the  computed  rms  values,  and  counting  the  number  of  runs  (crossings  above  and 
below  the  median  value)  in  the  set  of  rms  values.  The  number  of  runs  is  then 
compared  at  the  5-percent  level  of  significance  to  that  for  theoretical  Gaussian 
distribution,  and  the  hypothesis  that  the  data  are  stationary  is  either  accepted 
or  rejected.  Printouts  list  the  average  value,  maximum  absolute  value,  and  rms 
value  for  the  data  in  the  same  blocks  in  which  the  PSDs  will  be  computed  so 
inconsistent  blocks  of  data  can  also  be' deleted  from  the  analyses. 

The  statistical  moments,  or  weighted  excursions  about  the  mean,  were  also 
computed  as  part  of  the  verification  process.  Only  the  first  four  moments  are 
generally  of  interest  and  are  described  below: 

-  1/n  2x^  -  x 

-  l/nS-CXj^-x)2 

-  l/n'2.(x^-x)^ 

-  1/nKx^x)4 

The  first  two  moments  are  commonly  referred  to  as  the  mean  and  variance, 
while  the  third  and  fourth  moments  are  sometimes  called  the  skewness  and 
kurtosis  of  the  distribution.  More  often,  a  normalized  version  of  the  skewness 
and  kurtosis  are  computed  as  shown  below: 

Skewness  - 


Kurtosis  -  M, 

4 

*’**2 

<m2) 

A  normal  distribution  ha3  a  skewness  of  zero  (as  do  all  symmetric 
distributions)  and  a  kurtosis  of  3.  If  the  kurtosis  of  the  data  set  is  higher 
than  this,  it  is  an  indication  that  there  is  considerable  data  out  at  high 
amplitude  values.  Because  the  kurtosis  value  is  based  on  the  fourth  power  of 
the  distance  from  the  mean,  it  is  extremely  sensitive  to  wild  points,  even  in 
small  quantities.  This  value  has  proven  to  be  an  effective  screen  against 
abnormal  data. 
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APPENDIX  I.  VIBRATION  DATA  ANALYSIS 


The  acceleration  power  spectral  density  was  computed  by  dividing  the  time 
domain  data  into  successive  blocks  of  2048  points,  converting  to  the  frequency 
domain  using  the  Fast  Fourier  Transform  (FFT),  multiplying  this  result  by  its 
complex  conjugate,  and  linearly  accumulating  (averaging)  these  results  over  the 
entire  data  run.  The  number  of  linear  averages  applied  was  dependent  upon  the 
amount  of  valid  data  available  for  a  particular  data  run.  For  a  linearly 
averaged  process,  the  number  of  statistical  degrees-of-freedom  is  equal  to  twice 
the  number  of  averages  used.  The  amount  of  averaging  (degrees-of-freedom) 
determines  the  degree  of  confidence  that  the  value  measured  is  a  true 
representation  of  the  actual  physical  phenomena.  An  error  band,  based  on  the 
number  of  averages,  can  be  computed  for  various  confidence  levels  from  the 
chi  squared/degrees-of-freedom  distribution.  When  the  FFT  algorithm  is  used,  an 
assumption  is  made  that  the  time  record  being  transformed  (data  block)  is 
repeated  throughout  time.  If  the  time  record  contains  an  integer  number  of 
cycles  within  the  data  block,  the  assumption  is  valid  and  the  waveform  is  said 
to  be  periodic  within  the  time  record.  In  most  cases,  the  data  are  not  periodic 
within  the  record  which  causes  a  truncation  of  the  signal  at  the  end  of  the  data 
block.  Since  the  assumption  is  one  of  a  repeated  waveform,  the  analysis  process 
assumes  that  the  truncation  is  repeated  throughout  the  entire  data  record.  The 
effect  in  the  time  domain  is  an  apparent  discontinuity  in  the  representation  of 
the  data  signal.  In  the  frequency  domain,  the  discontinuity  appears  as  side 
lobes  or  additional  frequency  components  and  is  known  as  leakage.  A  time  domain 
truncation  technique  known  as  windowing  is  employed  to  reduce  the  leakage.  A 
Hanning  window  was  employed  for  this  analysis. 

The  dynamic  characteristics  of  a  physical  system  can  be  described  by  a 
frequency  response  function,  H(f) ,  which  relates  system  response  to  a  given 
Input  by: 

Y(f)  -  H(f)  X(f) 


where 


Y(f)  -  Fourier  transform  of  system  response 

H(f)  -  frequency  response  function 

x(f)  -  Fourier  transform  of  system  input. 

Since  both  Y(f).  and  X(f)  can  be  complex  values  (real  and  imaginary 
components),  H(f)  can  be  a  complex  value.  In  practical  terms,  the  relationship 
between  the  real  and  Imaginary  components  at  any  frequency  are  used  to  find  the 
phase  relationship  at  that  frequency. 


The  frequency  response  function  is  computed  as  the  ratio  of  the  cross  power 
spectrum,  G  ,  divided  by  the  input  auto  power  spectrum,  G  .  The  cross  power 
spectrum  is^lefined  as  the  Fourier  transform  of  Y(t),  Y(f),  multiplied  by  the 

complex  conjugate  of  the  Fourier  transform  of  x(t),  X*(f) . 

H(f)  -  ^(O/G^f) 

where 


H(f)  -  frequency  response  function 

G  (f)  -  cross  power  spectrum  between  input  and  response 
(£(f)  -  auto  power  spectrum  (PSD)  of  input. 

Since  the  frequency  response  function  contains  real  and  imaginary  values, 
it  can  be  difficult  to  use  in  normal  form.  For  this  analysis,  the  polar  form 
(magnitude  and  phase)  was  calculated  and  presented  (magnitude  only).  The  polar 
form  was  calculated  from: 

H(f )  -  <Hr(f)2  +  H^f)2)172 


where 


H(f)  -  magnitude  of  frequency  response  function 

H  (f)  -  real  part  of  frequency  function 

H^(f)  -  imaginary  part  of  frequency  response  function. 

When  calculating  the  frequency  response  function,  extraneous  inputs  and 
system  nonlinearities  can  cause  errors  in  the  computation.  The  coherence  / 

function  is  a  procedure  used  to  check  the  validity  of  the  frequency  response 
function  and  thus  can  be  used  to  check  causality  (input  a  causes  response  b) 
between  two  signals.  The  value  of  the  coherence  function  is  a  dimensionless 
number  between  zero  and  one.  A  zero  value  indicates  that  the  system  may  have 
extraneous  inputs  or  may  be  highly  nonlinear.  A  value  of  1  indicates  complete 
coherence  between  input  and  output  (linear  system,  only  one  input).  The 
coherence  function  was  calculated  as  follows: 

cVf>  -  vf>2/(c~<f>  vf» 

where 

2 

C  (f)  -  coherence  function 
G  X^f)  ”  magnitude  of  cross  power  spectrum 
G^rx(f)  -  input  auto  power  spectrum 
G^(f)  -  response  auto  power  spectrum. 

.'he  computational  techniques  used  to  produce  the  cross  power  spectrum, 
input  auto  power  spectrum  and  response  power  spectrum  (block  size,  averaging, 
windowing,  etc.)  were  the  same  as  for  the  P2D  computations  discussed  previously. 
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These  functions  were  not  computed  directly;  however,  the  program  dvpsd  was 
used  to  compute  the  transmissibility  function  described  here  as: 

T  (f)  -  transmissibility  function 
GX^(f)  —  input  auto  power  spectrum 
Gy/f >  ~  resPonse  auto  power  spectrum. 

It  can  be  shown  that  the  transmissibility  function  is  equal  to  the  square 
of  the  magnitude  of  the  frequency  response  function  if  the  coherence  function 
has  a  value  of  one  and  is  larger  than  that  in  all  other  cases.  Thus,  a 
pseudo- transfer  function  was  calculated  by  determining  the  transmissibility  and 
taking  the  square  root  of  the  result.  This  pseudo- transfer  function  represents 
an  upper  bound  of  the  actual  frequency  response  function. 
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APPENDIX  J.  DEVELOPMENT  OF  LABORATORY  VIBRATION  TEST  SCHEDULES 


A  series  of  computer  programs  have  been  written  for  the  purpose  of  creating 
laboratory  vibration  schedules  from  field  data.  The  programs  are  designed  to 
automate  the  process  as  much  as  possible  and  require  operator  intervention  only 
when  judgements  are  required  (such  as  shaping  the  schedules).  The  following 
description  is  a  "walk  through"  of  a  typical  development  project  once  the  field 
(raw)  data  have  been  collected  and  verified. 

The  program  PSDFL  computes  PSDs  from  raw  data  and  stores  them  in  a  file. 
The  program  assumes  that  the  data  are  in  normal  PCM  digital  multiplexed  form  and 
are  located  on  the  data  disc.  The  program  makes  use  of  a  control  file,  $PSDCL, 
which  contains : 

a.  The  name  of  the  calibration  file  (a  file  which  contains  the  values 
required  to  convert  the  integer  data  to  engineering  units). 

b.  The  name  of  the  amplitude  distribution  file  (a  file  which  contains 
information  on  the  location  of  the  desired  channels  in  the  data  stream  and  a 
brief  location  description  of  the  desired  channels). 

c.  The  data  sample  rate. 

d.  The  FIT  block  size  (currently  limited  to  2048). 

e.  The  word  length  (number  of  bits)  of  the  data  (typically  10). 

f.  The  lower  and  upper  frequencies  to  be  analyzed. 

g.  The  use  (or  not)  of  a  Hanning  window. 

When  the  program  is  run,  the  operator  enters  the  run  number,  the  starting 
track  (location  of  data  on  the  disc)  and  the  number  of  tracks  to  be  used 
amount  of  data  analyzed) .  The  program  creates  a  data  storage  file  with  a 
based  on  the  run  number.  Each  segment  of  time  domain  data  is  windowed  with 
a  Hanning  window  to  reduce  spectral  leakage  if  this  option  is  selected.  The 
program  then  computes  three  PSDs  for  each  channel  listed  in  the  amplitude 
distribution  file.  The  first  PSD  computed  is  a  linear  average  of  PSDs  computed 
for  segments  of  data  of  a  length  defined  in  the  control  file.  During  this 
analysis  period,  the.  standard  deviation  and  the  peak  value  for  each  spectral 
line  are  computed.  One  standard  deviation  Is  then  added  to  the  average  value  to 
create  a  PSD  which  represents  the  average  value  plus  one  standard  deviation  at 
each  spectral  line  (limited  by  the  peak  value  in  each  spectral  line) . 
Experience  has  shown  that  the  variation  in  the  data  over  a  test  run  is 
sufficient  to  cause  the  standard  deviation  to  be  slightly  larger  than  the 
average  value,  thus  adding  one  standard  deviation  to  the  average  has  the  effect 
of  doubling  the  average  value.  Both  this  PSD  (average  plus  one  standard 
deviation)  and  the  peak  value  PSD  (largest  value  at  each  spectral  line)  are  also 
saved  in  the  file. 
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Each  data  file  corresponds  to  a  test  run  (specific  test  course,  speed,  and 
test  condition)  and  consists  of  a  major  title  (a  line  that  specifies  the 
amplitude  distribution  file  used)  and  information  on  each  channel  analyzed, 
organized  as  follow: 

a.  An  ASCII  title  line  which  contains  the  run  number  and  a  brief 
description  of  the  channel  location. 

b.  An  ASCII  information  line  which  contains: 

(1)  The  FFT  block  size. 

(2)  The  data  sample  rate. 

(3)  The  analysis  bandwidth. 

(4)  The  peak  value  contained  within  the  PSD. 

(5)  The  number  of  linear  averages  used. 

(6)  The  use  (or  not)  of  a  Hanning  window. 

c.  A  binary  record  block  which  contains  formatted  integer  data.  (The 
data  are  formatted  by  dividing  the  entire  spectrum  by  the  peak  value,  then 
multiplying  these  values  by  10,000  to  save  file  storage  space;  i.e.,  integer 
numbers  use  1/2  the  storage  space  of  real  numbers.) 

The  program  PSDFL  is  run  on  the  data  for  N  runs  with  data  from  X  channels 
being  analyzed  from  each  run  (note  that  X  need  not  be  the  same  for  each  run) . 
The  program  will  produce  N  PSD  files  (one  for  each  run) . 

At  this  point,  an  overlay  process  is  started.  A  program  called  PKPSD  is 
used  to  create  the  overlayed  spectrum.  The  operator  enters  (via  a  control  file) 
the  name  of  the  master  storage  file,  the  names  of  the  data  files  created  by 
PSDFL  (essentially,  the  data  runs  to  be  analyzed) ,  the  data  elements  (channel 
locations)  within  those  files  to  be  analyzed  (those  locations  for  which  a  common 
spectrum  is  desired) ,  which  PSD  values  are  to  be  used  in  the  overlay  process 
(average,  average  +1  standard  deviation  or  peak)  and  the  number  of  standard 
deviations  (due  to  the  overlay  process)  to  be  added  to  the  average  value 
computed  during  the  overlay  process .  The  program  does  not  require  that  the  data 
elements  to  be  analyzed  oe  the  same  for  each  file  analyzed;  however,  all  data  to 
be  compared  must  have  the  same  analysis  bandwidth  (a  check  is  made  by  the 
program  and  the  data  are  rejected  if  this  is  not  the  case). 


The  data  normally  used  for  the  overlay  process  for  e»ch  particular  channel 
to  be  included  in  the  process  are  the  average  plus  one  standard  deviation  PSD 
data.  Some  of  the  rationale  for  the  use  of  these  data  (rather  than  the  average 
or  the  maximum)  are  as  follow: 

a.  The  paak  value  at  each  spectral  point  is  a  value  that  occurs  rarely 
(once)  and  is  susceptible  to  contamination  by  noise  (false  data  from  any 
source).  To  create  a  schedule  based  on  peak  data  is  to  assume  that  the  test 
item  is  subjected  to  such  an  environment  for  a  reasonable  length  of  time,  which 
is  false.  The  result  is  an  overtest. 

b.  Some  conservatism  is  built  in  by  adding  a  standard  deviation  to  the 
average  to  account  for  other  terrain  and  vehicles  for  which  no  data  were 
obtained,  effects  of  tire  pressure,  vehicle  wear,  track  tension,  etc. 

The  program  computes  the  average  and  the  standard  deviation  of  all  PSD's  to 
be  overlayed  (at  each  spectral  line)  and  again  saves  the  average  plus  X 
(selected  by  the  operator,  but  normally  1)  standard  deviations  (limited  by  the 
largest  value  at  each  spectral  line  encountered  during  the  overlay  process). 
Conservatism  is  thus  added  twice  during  the  overlay  process.  The  use  of  the 
average  plus  standard  deviation  curve  from  the  individual  location  PSD's 
accounts  for  variations  in  a  single  location  during  a  data  run  while  the 
addition  of  a  standard  deviation  during  the  overlay  process  accounts  for 
variations  among  locations  or  test  runs  (courses) .  The  standard  deviations  are 
computed  during  separate  processes  and  do  not  have  the  same  values. 

It  is  the  responsibility  of  the  operator  to  determine  the  data  locations 
and  data  runs  to  be  overlayed  and  to  insure  that  the  data  for  each  of  the  runs 
and  locations  chosen  are  valid.  The  program  PKPSD  prints  the  channel 
description  of  each  channel  used  in  the  overlay  process  so  that  the  operator  may 
verify  that  the  requested  data  have  been  analyzed.  The  rms  value  of  each  of 
the  individual  spectra  is  computed  and  printed  with  a  recommendation  to  discard 
any  data  with  an  rms  value  less  than  1/3  the  average  (all  elements  overlayed)  or 
more  than  3  times  the  average  rms  value.  Locations  which  fit  this  description 
(if  any)  are  Identified,  but  no  action  to  delete  these  locations  is  taken  by  the 
program. 

A  program  called  CRVIB  is  next  used  on  the  master  (overlayed)  spectrum  to 
define  the  break  points  (frequency  and  amplitude  points  which  define  the 
spectrum  in  the  digital  vibration  control  system)  and  allow  the  operator  to 
smooth  the  spectrum  for  laboratory  simulation.  The  master  spectrum  is  plotted 
on  a  graphics  terminal  in  a  log- log  format,  and  the  cursor  is  turned  on.  The 
operator  may  position  the  cursor  at  any  location  on  the  plot  and  enter  the 
location  of  the  point  (frequency  and  amplitude)  by  hitting  any  key  (except  E). 
The  location  of  the  cursor,  in  screen  coordinates,  is  converted  to  engineering 
units  of  frequency  (Hz)  and  amplitude  (g  /Hz) ,  and  these  values  are  stored  in  a 
file.  The  program  utilizes  the  "rubber  band  line"  technique  when  the  cursor  Is 
moved,  thus  giving  the  operator  a  preview  of  the  smoothed  spectrum.  The  use  of 
a  log- leg  domain  plot  allows  better  resolution  in  the  low  frequency  portion  of 
the  spectrum  which  contains  most  of  the  spectral  energy. 
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Laboratory  test  time  can  be  reduced  over  actual  exposure  time  in  the  field 
by  applying  an  exaggeration  factor  to  the  amplitude  (PSD)  values  of  the 
schedule.  The  rationale  is  based  on  an  equivalent  damage  theory  assuming  that 
the  failure  mechanism  is  fatigue.  The  techniques  for  computing  an  exaggeration 
factor  are  based  on  a  detailed  knowledge  of  the  test  vehicle  real-world  scenario 
(i.e.,  the  real-world  exposure  time),  the  desired  test  time  and  the  material 
fatigue  characteristics.  The  accepted  equation  for  test  time  compression  is: 

oy V 2)b/n  -  Tj/Tj. 

where 

-  field  PSD  amplitude 
Uj  -  laboratory  PSD  amplitude 
T^  -  real  world  exposure  time 
Tj  “  laboratory  test  time 
b  -  9  (endurance  curve  con. .ant) 
n  -  2.4  (stress -damping  constant). - 

The  breakpoints  (and  thus  the  spectrum)  are  now  exaggerated  using  the 
program  CRPSD.  It  is  left  to  the  operator  to  compute  the  value  of  the 
exaggeration  factor  based  on  the  desired  test  time  to  real  time  ratio.  The 
operator  enters  the  exaggeration  factor  (by  which  the  PSD  levels  arc  amplified) 
and  the  name  of  a  new  storage  file  into  which  the  exaggerated  spectral  data  will 
be  stored.  The  program  takes  the  breakpoint  data  from  the  applicable  file, 
multiplies  the  amplitudes  by  the  exaggeration  factor  and  stores  the  values  Into 
the  new  file.  The  program  then  takes  the  new  breakpoints  (exaggerated)  and 
generates  a  complete  spectrum  (with  the  same  analysis  bandwidth)  by  calculating 
points  between  breakpoints  using  a  straight  line  fit  in  the  log-log  domain. 
This  created  spectrum  (the  actual  laboratory  schedule)  is  stored  in  a  separate 
file  also  named  by  the  operator.  A  program  called  SHAKR  is  then  run  using  the 
just  developed  laboratory  spectrum  as  data  to  insure  that  shaker  displacement 
limits  are  not  exceeded  by  the  test  schedule.  The  program  compares  the  spectral 
value  at  each  frequency  from  5  to  15  Hz  with  the  maximum  value  known  to  function 
properly  on  the  USACSTA  system  and  alerts  the  operator  to  any  changes  required 
In  the  spectrum. 
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